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We present a theory of the elementary spin excitations in transition metal ferromagnet nanopar- 
ticles which achieves a unified and consistent quantum description of both collective and quasipar- 
ticle physics. The theory starts by recognizing the essential role played by spin-orbit interactions in 
determining the energies of ferromagnetic resonances in the collective excitation spectrum and the 
strength of their coupling to low-energy particle-hole excitations. We argue that a crossover between 
Landau-damped ferromagnetic resonance and pure-state collective magnetic excitations occurs as 
the number of atoms in typical transition metal ferromagnet nanoparticles drops below approxi- 
mately 10^, approximately where the single-particle level spacing, 5, becomes larger than , y/aErcs, 
where E-^cs is the ferromagnetic resonance frequency and a is the Gilbert damping parameter. We 
illustrate our ideas by studying the properties of semi-realistic model Hamiltonians, which we solve 
numerically for nanoparticles containing several hundred atoms. For small nanoparticles, we find 
one isolated ferromagnetic resonance collective mode below the lowest particle-hole excitation en- 
ergy, at ErcB ~ 0.1 meV. The spectral weight of this pure excitation nearly exhausts the transverse 
dynamical susceptibility spectral weight. As 5 approaches ^/aEias, the ferromagnetic collective ex- 
citation is more likely to couple strongly with discrete particle-hole excitations. In this regime the 
distinction between the two types of excitations blurs. We discuss the significance of this picture 
for the interpretation of recent single-electron tunneling experiments. 



I. INTRODUCTION 



In bulk condensed matter systems, normal metals are 
Fermi liquids; their low-energy excitation spectra can be 
placed in one-to-one correspondence with those of corre- 
sponding non-interacting electron systems as argued by 
Landau more than 50 years ago. Recent single-electron 
tunneling spectroscopy studies of metallic nanoparticles 
, in which the discrete excitation spectra of small sys- 
tems containing less than one thousand to tens of thou- 
sands of atoms are investigated, have allowed Landau's 
enormous simplification of interacting fermion physics to 
be examined quite directly. With a few caveats and some 
exceptions, the wide variety of interesting phenomena 
that have been studied using single-electron tunneling 
spectroscopy can be understood using ideas from inde- 
pendent particle quantum mechanics. Although they can 
still be regarded as Fermi liquids for many purposes, met- 
als with nearly continuous broken symmetries (in particu- 
lar the ferromagnetic transition metals that are the focus 
of this article) support low-energy collective excitations 
in addition to Landau's particle-hole excitations. When 
spin-orbit coupling is neglected, the broken symmetry 
of itinerant electron ferromagnets is continuous and the 
collective excitations are Goldstone bosons whose energy 
vanishes in the long-wavelength limit. Recent single- 
electron tunneling spectroscopy studies 0,0 have for the 
first time succeeded in resolving the excitation spectra 
of ferromagnetic transition metal nanoparticles with di- 



ameters below 4 nm. The discrete resonances seen in 
the tunneling experiments measure the low-energy many- 
electron excitation spectra of a single-domain ferromag- 
netic nanoparticle. The ultimate objective of this paper 
is to shed light on the physics that is responsible for the 
rich and complex behavior seen experimentally, which 
includes hysteretic behavior, non-monotonic field depen- 
dences, and a much denser low-energy excitation spec- 
trum than would be expected based on a non-interacting 
quasiparticle model. Our approach is based on a weak- 
coupling description of a metallic ferromagnet in which 
spin-orbit interactions cause collective and particle-hole 
excitation to be coupled at low-energies, and the classi- 
cal micromagnetic description of a ferromagnetic metal 
nanoparticle emerges naturally when quantum effects are 
unimportant. 

Our theory builds on earlier work 0, which cap- 
tures several features of the experimental spectra, espe- 
cially when non-equilibrium excitationsQ are considered, 
but does not provide a unified and consistent quantum 
description of how collective and quasi-particle excita- 
tions are coupled by spin-orbit interactions. The purpose 
of this paper is to develop such a description and illus- 
trate its implications by applying a simplified but quali- 
tatively realistic microscopic model that we have recently 
introduced0- 

The elementary spin excitations in bulk itinerant- 
electron ferromagnets are of two kinds: collective spin 
excitations (spin-waves) and spin-flip particle-hole exci- 
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tations (Stoner excitations). Spin- wave excitations are 
related to the collective magnetization degree of freedom 
and form a branch m q — lu space, which is gapless if the 
system is isotropic in accord with Goldstone's theorem. 
The main effects of spin-orbit interactions in the bulk 
are to generate an energy gap E^^s in the q = collective 
mode, which is of the order of the consequent magneto- 
crystalline anisotropy energy per atom - K ^ 0.1 meV 
in Cobalt - and to introduce the possibility of decay 
of long-wavelength collective excitations into particle- 
hole excitationsj^, 0|, a process that contributes sub- 
stantially to the collective excitation lifetime [s^l for the 
case of NiFe thin films. In the absence of spin-orbit in- 
teractions, the ferromagnetic resonance is coupled only 
to spin-flip particle-hole excitations which at long wave- 
lengths have a gap of the order of the spin-splitting field 
A. Gapless spin- flip particle- hole excitations are possible 
only at wavevectors exceeding the minimum q-space sep- 
aration between majority and minority spin Fermi sur- 
faces. The separation in energy at long-wavelengths be- 
tween collective modes and spin-flip particle-hole contin- 
uum implies that low-energy collective modes are only 
weakly damped. Beyond a critical value of q the spin- 
wave branch merges with the continuum. Thus spin- 
waves can decay into Stoner excitations. The strength of 
this decay process is sensitive to the character of the or- 
bitals involved in the particle-hole excitations llj. In the 
present work we will focus on how this description of ele- 
mentary spin excitations in itinerant ferromagnets has to 
be altered when the level spacing for quasi-particle exci- 
tations 6 is finite and approaches i^res, a condition that is 
satisfied in the nanometer particle-size range. Note that 
S is inversely proportional to particle volume, while -Ejcs 
is approximately volume independent for large nanopar- 
ticles. Since in a finite system there is no wavevector 
conservation, collective modes and spin-flip particle-hole 
pairs cannot be simply separated, unlike the bulk case. 
We will show that this fact, together with the essential 
role played by spin-orbit interaction, has profound conse- 
quences on the nature of the elementary spin excitations 
in ferromagnetic metal nanoparticles. 

The paper is organized in the following way. In SecHTI 
we introduce two similar microscopic models for a mag- 
netic metal nanoparticle and explain how they are re- 
lated to the phenomcnological model consider previously. 
The models differ in that one accounts for the difference 
in the strength of exchange interactions between s and 
d electrons in transition metals, a feature whose conse- 
quences we wish to address specifically. In Sec. IIIII we 
derive a path integral formulation of theories based on 
these models This point of view provides a conve- 
nient language for explaining the interplay between col- 
lective modes and particle-hole excitations, and for mak- 
ing contact with classical micromagnetic theory. The 
spin-orientation fluctuation propagator in the Gaussian 
approximation is discussed in Sec. IIVI The poles of the 
fluctuation propagator occur at the elementary spin ex- 
citations of the system. We will show that the ferro- 
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FIG. 1: Crossing of relevant energy scales as a function of the 
number of atoms J\fa in a magnetic nanoparticle. Here 5 is 
the single-particle mean level spacing; _E,os is the energy of 
the coherent spin collective mode or ferromagnetic resonance 
energy; K is the magneto-crystalline anisotropy energy per 
atom, and a is the Gilbert damping factor. 



magnetic resonance energy iJres can be expressed as the 
quotient of anisotropy energy and Berry curvature coeffi- 
cients which specify the Gaussian expansion of the action 
at low frequencies. We then discuss how the resonance 
evolves with particle-size. Our main conclusions is sum- 
marized in Fig. ^ For large particle-sizes the ferromag- 
netic resonance weight is distributed over a large number 
of particle-hole excitations, while for smaller particle sizes 
the ferromagnetic resonance appears as a pure quantum 
excitation. The crossover between the two regimes oc- 
curs approximately where the level spacing S is equal to 
\/a Ei-Qs where a is the bulk resonance's Gilbert damp- 
ing factor[l3j. For typical transition metal nanoparticles 
this condition is satisfied in particles containing of order 
10^ atoms. For smaller particles, avoided crossings be- 
tween collective and individual particle-hole excitations 
will occur occasionally as a function of system parame- 
ters, for example as a function of an external magnetic 
field used to reorient the magnetization. For nanoparti- 
cles containing fewer than of order 10^ atoms, avoided 
crossings with particle-hole excitations will usually not 
occur at any field orientation, and the nanoparticle can 
be considered to be a molecular magnet in which only 
spin-orientation degrees of freedom are important at low 
energies In Sec. we will present numerical results 
for a few-hundred atom nanoparticle, which illustrate 
some of the these points. Finally in Sec. IVII we sum- 
marize our findings and comment on their relevance in 
understanding current tunneling experiments. 
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II. QUANTUM MODELS OF A 
FERROMAGNETIC METAL NANOPARTICLE 

In this article we consider two slightly different quan- 
tum models chosen to describe both collective and quasi- 
particle physics in a magnetic metal nanoparticle. We 
will denote them as Local D-orbital Exchange model 
(LDE model) and Long-Range Exchange model (LRE 
model) respectively, for reasons that will become clear be- 
low. We will see that these models, when solved within a 
mean-field approximation, are essentially equivalent and 
provide a convenient quantum description of a ferromag- 
netic nanoparticle when the magnetization is coherent 
(spatially constant across the sample). Our use of these 
models is motivated partly by the evident success of spin- 
density-functional theory in describing ferromagnetism in 
transition metals; our formalism could easily be adapted 
to be compatible with this method of calculating the en- 
ergies of different magnetic configurations. The models 
are intended to be sufficiently realistic to capture generic 
aspects of transition metal nanoparticle magnetism, but 
evidently miss features that can be important in practice 
such as variation in exchange interaction strength and 
inter-atom hopping amplitudes near the surface of the 
nanoparticle. 



is written in terms of creation and annihilation op- 
erators c] ,, „ and c. ,, „ labeled by atomic-site indexes 
atomic angular momentum indexes fJ.i,H2 and spin 
indexes s,s'. We choose the spin-quantization axis to 
be along the direction of the magnetization, which is 
specified by a unit vector Cl{Q, $) where Q and <& are 
the usual angular coordinates defined with respect to 
the f.c.c. crystal axes. The parameters i^^^^^^s are 
Slater-Koster parameters obtained after performing 
a Ldwdin symmetric orthogonalization procedure 16] on 
the set of Slater-Koster parameters for non-orthogonal 
atomic orbitals of bulk spin-unpolarized Co[l3|. Here the 
exchange term is a short-range spin-interaction involving 
only the electrons spins of d-orbitals on the same atomic 
site: 

^^cxch = —'2'Udd 2, ^d,i ■ Sd,i , (3) 



where 



A. Local D-orbital Exchange Model 

The first model that we consider accounts qualita- 
tively for the orbital dependence of exchange interaction 
strengths in a transition metal itinerant-electron ferro- 
magnets 7|. We model the nanoparticle as a cluster of 
Afa atoms located on the sites of a truncated crystal. The 
numerical results we present here are for a cobalt cluster 
whose truncated f.c.c. crystal is circumscribed by a hemi- 
sphere whose equator lies in the XF-plane of the f.c.c. 
crystaljU- The choice of a hemisphere is motivated by 
the tunneling experiments of Ref. |3, 0- We use a s-p-d 
tight-binding model for the quasiparticle orbitals, with 
18 orbitals per atom, including the spin-degree of free- 
dom. Nine orbitals per Co atom are occupied in neutral 
nanoparticles. The full second-quantized Hamiltonian is, 

W = ^band + -ffcxch + ^^so + Hzcc ■ (1) 

The second-quantized one-body term 

i,j s fj,i,fj,2 

I 



a being a vector whose components cr", a = x,y,z 
are the three Pauli matrices. The parameter Udd in 
Eq. determines the strength of the exchange inter- 
action and we will set it equal to 1 eV in our numerical 
calculations 7] . This choice leads to the correct magneti- 
zation per atom in the bulk. 

The spin-orbit coupling Hgo is a local one-body opera- 
tor 

where the atomic matrix elements (/i, s|L • S\^',s') = 
{i, fj,, s\L- S\i, fi' , s') depend on the spin-quantization axis 
specified by the angles 8 and $. The energy scale £,d, 
which characterizes the coupling between spin and or- 
bital degrees of freedom, varies in the range from 50 to 
100 meV in bulk 3d transition metal ferromagnets[l8|. 
Finally, Hzcc is a local one-body operator describing the 
Zeeman coupling of the orbital and spin degrees of free- 
dom to an external magnetic field Hcxt'- 
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with gs = 2. Hzcc plays an important role in manipulat- 
ing many-body states. In Ref. (0) we have investigated 
the spectrum of the microscopic Hamiltonian of Eq. 
treating the quartic exchange interaction in the mean- 
field approximation. We have focused in particular on 
the mesoscopic physics of the quasiparticle energies and 
their complex behavior as a function of the magnetiza- 
tion orientation and external magnetic field orientation 
and strength. This analysis, although very relevant to the 
understanding of tunneling experiments, does not tell us 
anything about the quantization of the collective magne- 
tization orientation dynamics. For this purpose and for 
the purpose of making a connection with the classical mi- 
cromagnetic theory, the path integral approach described 



in the next section provides a more useful language. 



B. Long- Range Exchange Model 

Our microscopic LDE model, when solved in the mean- 
field approximation, is related to a toy model, originally 
introduced in Ref. which can be regarded as the sim- 
plest possible model of a ferromagnetic metal nanopar- 
ticle. The toy model Hamiltonian assumes identical ex- 
change constants between all pairs of single-particle or- 
bitals in the nanoparticle: 



n,m s,s' ,t,t' 



S . 



(7) 



In Eq. 13 cjj g and c„ ^ are Fermion creation and anni- 
hilation operators for a quasi-particle state character- 
ized by orbital energy e„ and spin component s; S — 
\^n'^n s'^s' ,sCn^s 1^ the total Spin of the nanoparticle. 
The single particle orbitals will have an average spacing 
inversely proportional to the volume of the nanoparticle 
(or the number of atoms Ma) and are expected to exhibit 
spectral rigidity The one-body term in this Hamilto- 
nian, -ffband, should bc thought of as including a mean- 
field approximation to those spin-independent portions 
of the interaction not captured by the exchange term 
Hcxch = ■ S . The many-particle spectrum of this 

Hamiltonian has been discussed in detail in Ref. Q). 

This toy model Hamiltonian can be further 
augmented [l2| by a one-body spin-orbit coupling 
term Hso- We can write Hso asQ 

n,m,s 

with 

7,'* = iv^ Y = -v" v" =0 (9) 

where the conditions on the matrix elements ^ speci- 
fied in Eq. ^ ensure that Hgo is hermitian and invariant 
under time-reversal. 

Consider the orbital part of the microscopic Hamilto- 
nian given in Eq. (Q. Since i?band is quadratic, it can be 
diagonalized by a canonical transformation: 

-ffband — ^fli,fj.2,s'^i,fii,s'^j,fj.2,s ~ '^li.s'^n.s^n , 

S /^l,/12 71. S 

(10) 



where 

and |n, s) are the orthonormal eigenvectors. The eigen- 
values e„ in Eq. ifTUIl can be identified with the doubly 
degenerate orbital energies of the toy model Hamiltonian 
in Eq. iQ. If we now use Eq. Hll|) and the complete- 
ness of the eigenstates {|«, /i, we can rewrite the 
exchange interaction of Eq. Q in the form 

H....-~l^J-S--^}T.S^-S^^ (12) 

with Si = J^fi'^ifj.- From Eq. (|12|l it is clear that the 
toy model Hamiltonian with equal exchange constants 
is equivalent to a microscopic model with a long-range 
exchange interaction, coupling the spin of each atomic 
orbital at each atomic site to all others; hence the name 
of Long- Range Exchange model (LRE model). We em- 
phasize that the electron spins of all orbitals are cou- 
pled in this model, not just the spins of d-orbitals. From 
Eq. H12|) it follows that in a constrained mean-field ap- 
proximation where the local moment (5*^) is forced to be 
coherent across the sample and the spins of all atomic or- 
bitals are exchange-coupled, the microscopic LDE model 
and the the LRE model would be equivalent equivalent, 
(with 2Udd U/8). Differences arise when only the the 
d-orbitals in the LDE model are exchange-coupled, which 
we comment on below. 

Neither model includes any magnetostatic dipole- 
dipole interactions, which can be important in some cir- 
cumstances, for example when the nanoparticle is not 
close to spherical, but are easily incorporated in our dis- 
cussion. Note however, that both models lead to strong 
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shape dependence, because of surface effects which are 
important when the particle size is small. In Sec. (V) we 
will see that the structure of particle-hole excitations of 
the LDE model when only the spins of d-orbitals are cou- 
pled is richer than when quasiparticle majority and mi- 
nority spins are simply shifted by a rigid exchange field, 
as in the case of the LRE model. 



III. AUXILIARY FIELD FORMULATION 

A. Coherent state functional integral and 
Hubbard-Stratonovich transformation 

We consider first our LRE toy model. The extension 
to the LDE model is straightforward and we will com- 
ment on it below. Following some familiar steps poll2ll|. 



we write the interacting fermion partition function as an 
imaginary time coherent state path integral 



I?[^(t) V^t)] exp(-5) , 



(13) 



where the action S is 



dr 



n,s 

(14) 

Here /3 = l/k-gT, /i is the chemical potential and we use 
units such that h = 1. Since the exchange interaction 
term in our toy model is quadratic in the total electron 
spin - see Eq. Q - its contribution to the action for 
each time step Tk = ke = kp/N can be represented by a 
Gaussian integral over a real vector field = A(rfe) 



(15) 



Using this transformation for each time step, we obtain the following functional integral over a single auxiliary field 
A(r) that fluctuates in imaginary time 



I V[i^{r) ^{r)] I P[A(r)] exp(-5) , 



where the action S is 



S 



dr 



2U' 



A(r) • A(r) + ^ ^„^,(r)(— - fi)^PnA^) + i/ba„d(^(T), V^(r)) 



-Hso{tp{T),^p{T)) - -(A(t) +gsA*B-ffoxt) • ^ V'n,s(T)(7s,5'Vn,s(T) 



(16) 



(17) 



This action can be written in another form which is especially useful for small nanoparticles at low temperatures. Let 
us consider the fluctuating one-body Hamiltonian 



Hib{A{T)) = iJband + Hso - |(A(r) + gsfJ.BHcxt) ■ y^^cl.s'^s'.sC^,, 



(18) 



In Eq. Hlt)|l we can replace the functional integral over the fermionic coherent states iP{t),iIi{t) by a trace of the 
operator 



exp ( - / dTHu{Mr)) = exp ( - e^i?ib(Afe)) = H^^P ( - e-^ib(^fe)) ■ 

•'O 7. J, 



(19) 



By inserting at each time step the resolution of the identity / — |5'a(Afc))(^'a(Afc)|, where {\^ a{A.k))}a 
is a complete set of eigenstates of the one-body Hamiltonian J?ib(Ai.), and taking the T ^ limit, we obtain 
Z = / I?[A(r)] exp(-5), with 

S = ^Afe • Afc - In [(^'o(Ajv)| exp ( - eHib{KN-i)) |*o(AAr_i)) 



X . . . X (^-0(^1)1 exp ( - eHu{AN) I^-oIAat)) 



(20) 
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Note that the fluctuating field A(r) is bosonic and satisfies periodic boundary conditions A(0) = A(/3). Since e — > 
we can rewrite the action in the form 



S = 



I ^^(^^ ■ ^^^^ " 1 + e))|e-^-(^M|*o(A(r))) - 1 



dT 



^A(r) . A(r) + E,,{A{r)) + (vI/„|9*o/5A) • ^ 



(21) 
(22) 



Here Eu ( A(t)) is the ground-state energy of Hn, (A(t)) . 
The approximations that are implicit in these manipula- 
tions are valid when collective fiuctuations of the spin- 
splitting field are not strongly coupled to particle-hole 
excitations. 

Eq. 1)22(1 will play a crucial role in this paper. The 
quantity: 



Etot [A(7 



2U 



A(r).A(r)+i?ifc(A(r)) , (23) 



gives the quantum energy functional of the magnetic 
nanoparticle as a function of the fluctuating spin-splitting 
field A(t); its classical limit, obtained by evaluating 
at zero frequency (i.e. static A independent of time) 
Etot [A(t)] , corresponds to the phenomenological micro- 
magnetic energy functional for a coherent magnetic par- 
ticle with magnetostatic contributions neglected. 
The last term in the action 



5, 



Berry 



dT (*o|c'*o/aA) 



dAjr) 
dT 



(24) 



is a Berry phase term, which is related to the reduc- 
tion of the total spin component along the magnetiza- 
tion axis due to spi n deviations from the ground state 
configuration 122, E^- As we will see in Sec. l(IVp . the 
Berry phase contribution to the action captures the quan- 
tization condition of the collective elementary excitations 
and the way in which spin-orbit interactions affect this 
quantization condition. 

The Hubbard-Stratonovich decoupling can also be car- 
ried out in the case of the LDE model of Eq. O Here we 
should in principle introduce an auxiliary field Ad^i{T) 
that is dependent on the atomic site i. Complicated in- 
homogeneous and noncoUinear spin-splitting fields can 
indeed occur for very small nanoparticles with a few tens 
of atoms |23.l25| . In this paper however, we will consider 
only a coherent i.e. site-independent spin-splitting field, 
which is a good approximation for the relatively large 
nanoparticles {Afa > 50) that we are interested in. 



a time- independent spin-splitting field Amfj the Berry 
phase contribution to the action does not enter at this 
level. In the coherent-field approximation, the mean-field 
spin-splitting satisfies: 



, dEuiA) 
—Amf H — 

U OA A^ 

= —Amf — (^j\/F|<S'tot|^'AfF) 



= 0, 



(25) 



where |5'mf) = I^o(Amf))- Amf rnay be determined 
either by minimizing the energy functional or by solving 
the self-consistent-field equations implied by the second 
form of Eq. ( I25|l . The same set of mean-field equations 
can be derived directly from the more general expres- 
sion for the action in which the functional integral over 
fermion Grassmann variables is still present [2ll | . 



IV. GAUSSIAN FLUCTUATIONS 

Our theory of elementary magnetic excitations of 
metallic nanoparticles is based on a Gaussian fluctuation 
theory in which the action is expanded to second order 
around the static mean-field values of the auxiliary fields. 
It will prove informative to contrast two approaches that 
can be used to evaluate the Gaussian fluctuation ac- 
tion. In the first approximation we work directly with 
the Berry phase and energy terms in the action, which 
are rather transparently related to micromagnetic theory 
and Landau-Liftshitz semi-classical dynamics. This ap- 
proach cannot, however, deal directly with the coupling 
of collective and particle-hole elementary excitations. In 
the second approach we expand the action to second or- 
der around Amf, without introducing the quasi-static 
energy functional as an intermediate step. We will see 
that the two approaches are equivalent in the limit of en- 
ergies smaller than the minimum particle-hole excitation 
energy. 



B. Mean-Field Approximation 

The mean-field approximation for the action is ob- 
tained by finding the value of the spin-splitting field at 
which it is minimized. Since the minimum occurs for 



A. Energy Functional Approach 

If we are interested only in low-energy excitations, 
which have slow dependence on imaginary time, we can 
approximate the energy function by its dependence on 
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static field variations: 



imaginary time variable 



u 



A. 



A 







(26) 

where A" = A" - A^p, and E'mf = £^tot(AMF)- Since 
amplitude variations are energetically more costly, the 
dominant fluctuations of the order parameter will be ro- 
tations. The second derivative represents the expansion 
of the micromagnetic energy around its extremum; the 
dependence of energy on orientation is generally referred 
to as the anisotropy energy of the nanoparticle. 

The bosonic field A(t) is periodic in [0,(3]. There- 
fore the Berry phase contribution to the action given in 
Eq. (|24|l can be rewritten as the closed line integral that 
is path-dependent but not dynamics dependent: 



(27) 



C(A 



MP 



, . . . 

■/ dr[A^9^A2' - A^a^A^] , (31) 

"'0 



^Bciry 



which can also be written in frequency space as 

^Bcrry = X! C(Amf) [A^(-iw„) A^(lW„) 

-A^(~ZC^„) A^ilLUn)] . (32) 

Eq. 1)32(1 shows that the Berry phase contributes a term 
linear in frequency ujn to the imaginary part of the 
quadratic Gaussian fluctuation action. We will see that 
in the next section that the same contribution to the 
action arises from a small frequency expansion of the 
spin-fluctuation propagator kernel, defined in the next 
section. 



Since 

d 

(*o|-Y*o) 
oA 



oA oA 



dA 



r^'o|*o) 



_d_ 

'dK' 



;*oi*o> -0 



(28) 



the Berry phase iSBcny is pure imaginary. For small am- 
plitude rotations around the direction of Amf, the closed 
line integral Eq. (|27|) is equal to the product of the area 
of the enclosed path and the Berry curvature. 



C(A 



MP) 



z-Vx(*o 



dA' 



,9*0,9*0, ,9*0,9*0, 



'9A^'9A2'' 'dAv'dA^' 
(29) 

Unlike the integrand of Eq.( I27|l . the Berry curvature 
is gauge independent, i.e. it is independent of the ar- 
bitrary phase choice made for the many-electron wave- 
function at each magnetization orientation. For small 
fluctuations the integrand in Eq. H27|l can be chosen to 
be anything whose curl is the constant C(Amf)^. For ex- 
ample one convenient choice, analogous to the symmetric 
gauge choice for a constant magnetic field, leads to 



^Bcrry 



C(A 



MF 



■ (j) [z X A] ■ dA 



(30) 



if we choose the z direction to be the direction of Amf- 
This line integral can be parametrized in terms of the 



1. Calculation of the Berry's curvature 

We now briefl y d iscuss the numerical evaluation of 
Berry curvatures |22| evaluated at the mean-field-solution 
point. We start from the expression for 5Boriy given 
by Eq. lf?7jl . and consider paths on the unit sphere of 
exchange-field orientations. A small area closed path cen- 
tered on the mean-field orientation encloses a small area 
on the unit sphere whose normal is the direction of Amf, 
which we take as the z direction. Using Stoke's theo- 
rem and taking the Berry curvature out of the surface 
integral, we obtain 



^Bcrry — C(Amf) / dS^ ■ 
J Area 



(33) 



For our calculation we consider as a closed path a small 
right triangle in the xy-plane of sides A^^ = Ay = Amf^*, 
where is a small angle describing the rotation around 
the X and y axis. Then 



^Bcrry — -C(Amf)Amf^^ • 

This Berry phase can also be expressed asp2 



(34) 



i^Bcrry = -Im lu [(*MF | *0 ( A:, )) (*0 ( A^, ) | *o ( Ay )) (*o ( Ay ) | *mf) 



(35) 



where \'^o{Ax)) and |*o(Ay)) are the single Slater deter- minants of lowest energy for exchange fields that differ 
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from the mean- field value by Aa;2;, and Aj,y respectively. 
Eq. (|34|l can be obtained by expanding |4'o(A2:)) and 
|\E'o(Ay)) to quadratic order in Aj; and A^, remembering 
that (^'ol^^o) is pure imaginary. The expression for 
the Berry phase given in Eq. (|34|l is very suitable for nu- 
merical calculations, since the wavefunctions \'^o{Ax)) 
and (Ay)) can be easily calculated. The arbitrary 
phases of the wave functions at the three vertices ex- 
plicitly cancel since each wave function and its complex 
conjugate appears in Eq. If we choose the arbi- 

trary phases in such a way that the matrix elements of 
|\E'o(Aa;)) and |5'o(Ay)) with I^'mf) are real and positive, 
then the Berry phase is given by the simplified expression 



We obtain 



i^Borry = Im ln(«'o ( A^ ) | ^-Q ( Aj, )) 



(36) 



Finally, one can show that the Berry phase (or Berry 
curvature) of the many-body wave function is given by 
the sum of the Berry phases (or curvatures) of the occu- 
pied single particle states, even in the presence of spin- 
orbit interaction. The easiest way to prove this is to show 
that 

occ 

{^0\-^^0)^Y.(^n\-^<t>n), (37) 

oA „ oA 

where (/)„ are single-particle eigenstates of the one-body 
Hamiltonian. The gradient in Eq. (|37|l acts separately 
on the single-particle wave functions and hence it can be 
regarded as a sum of single-particle operators. It follows 
that the Berry phase can expressed as 



iSBerry = ^ Im ln(^„ (A^,) |(/i„ (Ay )) . (38) 



2. Berry phase without spin-orbit coupling 

For the case of no spin-orbit coupling, the Berry 
phase term can be calculated analytically and is equal 
to the number of singly-occupied orbitals times the usual 
spin-1/2 Berry phase. To show this, let us consider 
a many-body wave-function within the long-range ex- 
change model describing a state polarized in the direction 
n = A/A 



|vE'o(A)=n("(^)^sT+«(^)cIi) 11' 



it^iio) 



l*o(fi) , 
(39) 

where the index s runs over the — singly-occupied 
states and the index d over the 7V| doubly-occupied 
states. The functions u{Q) and v(fl) are written in terms 
of the angles 6 and <f> specifying the unit vector fi(8, $) 



u{n) = cos(e/2), 

v{h)=e"^ sin(e/2). 



(40) 
(41) 



n^T^ilO). (42) 



Hence, for rotations of the spin-splitting field 



1 - cos(e) - 

2sin(e) 



E 



- iV| 1 - cos(e) 

2 sin(e) 



$(43) 



and the Berry phase 



dr 



^d$(l-cos(e$)) 



dn 



(44) 



where A[f2] is the area enclosed by the path il(r) on 
the unit sphere, and A[r2]/2 is the usual spin-1/2 Berry 
phase. In the presence of spin-orbit interactions, the 
Berry phase becomes highly non-trivial as we discuss in 
SecEl 

For a path enclosing the small right triangle described 
in Sec. liV All A[n] = 9^/2 and 



5, 



Berry 



C(Amf) 



0^S 



' A^ 



(45) 
(46) 



where 5* = {N-\^ — Ni)/2 is the total spin in the nanopar- 
ticle ground state. 



B. Gaussian Fluctuations: Perturbation Theory 
Approach 

The Gaussian fluctuations theory can be derived more 
generally by formally integrating out the fermions and 
then expanding second order around Amf without mak- 
ing any quasi-static approximations. We discuss the case 
of the long-range exchange interaction first and then in- 
dicate what changes occur in the short-range (d-only) 
exchange model. The formally exact expression for the 
action is 

r^A-A-lndet(9, + i7ifc(AV + A")). (47) 

where Hn,{A) includes, in addition to its single-particle 
hopping and spin-orbit terms, the spin splitting terms 
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— (Amf + ^^oxt) • s + A"s". It is the second of these two 
terms that is treated perturbatively. First order terms in 
the expansion vanish, since the mean-field value of A is 
an extremum of the action. The second order terms can 
be obtained by a standard calculation with the following 
result: 

5fiuc = ^^A^(-iw„)i^^„(iw„)A"(ic^„) , (48) 

where the kernel a (*w„) is the inverse of the exchange- 
field-fluctuation propagator and is given by 

x{J\s''\I){I\s-\J) (49) 

where |/) is a mean- field electron eigenstate and ^/ is the 
corresponding eigenvalue with energies measured from 
the chemical potential. 

Fluctuations in the magnetization orientation around 
the mean field direction z = Amf/ Amf, are described 
by the transverse diagonal Kxx, Kyy and off-diagonal 
Kxy, Kyx components of the kernel. It is customary 
to introduce - (iun) and if_+(ia;„). When the spin- 
orbit interaction is present, if+ + (jw„) and K {itOn) do 

not vanish and play a role. Note that 



K+ + (-iw„) = /^+ + (iw„) . 



_(iw„) = K+ + {-~iuinY 
It follows form Eq. gSJl that 



(50) 



(51) 



(52) 



Kx X = Ma/U +{K+^+ +)/4 + (i<r+ + + _)/4 , 

(53) 

Kyy = Na/U +{K+^+ +)/4 -{K+++ _)/4 , 

(54) 



X y - i{K+ _ - if _ +)/4 - t{K+ + - if_ _)/4 , (55) 



K 
and 



Kyx = -i{K+--K-+)/A-i{K++-K--)l^. (56) 

The transverse fiuctuation action of the long-range ex- 
change model is quite simple in the absence of spin-orbit 

interactions. Firstly, the components if++ = K = 0, 

since the spin of the quasi-particle states is a good quan- 
tum number. (This property depends only on spin- 
rotational invariance and holds for the short-range ex- 
change model as well). Furthermore, since the the elec- 
tron mean-field eigenstates factorize into spin and orbital 



factors identical spin-up and spin-down wavefunctions 
are split energetically by Amf and we obtain 



(Amf + ^^cxt) 



lUJr, 



(57) 

In examining the consequences of this simple property for 
the Gaussian action kernel, we first consider Kxy{iuJn). 
We see from Eq. (|57|l that when there is no spin-orbit 
interaction, Kxy{ioJn = 0) =0, a property that reflects 
invariance under rotations of exchange-fleld orientation 
around the z axis. By expanding Kxy{ii-Un) around iujn = 
0, we obtain 



^^MF 



IS 



A2 
^MF 



(58) 



By inserting Eq. (|58|l into Eq. H48() and comparing 
the result with Eq. H32() . we can identify the term of 
KxyiiijJn) linear in iuJn, zS/A^p, with the Berry cur- 
vature, C(Amf), in agreement with Eq. 14t)|) . 

Similarly, the expansion of Kxxii^n) — Kyy{iujn) 
around iuJn — yields 



Kxx{ii^n) = Kyy{iu}n) ~ -JJ- 



2A 



MF 



= + O(c.^), 



(59) 



where in the last equality of Eq. we have used the 
mean-field relationship ([/ /Ma) S = Amf between the 
total magnetization of the nanoparticle and the spin- 
splitting field. The vanishing of the diagonal compo- 
nents Kxx and Kyy at = is again a result of ro- 
tational invariance. In the absence of spin-orbit interac- 
tions there is a collective mode at zero energy because of 
magnetization-orientation rotational invariance. 

The physics that we wish to investigate in this pa- 
per is largely contained in the way that the simple re- 
sults outlined above are altered by spin-orbit interac- 
tions. As we have stated, both Berry-phase and the 
energy-term in the action become highly non-trivial, a 
property that we now examine from the perturbation 
theory point-of-view. From Eqs. (|5U|I - H5()|I it is clear 
that Kxx{i^n), Kxx{uJn), Kxy{uJn) and Kxy{Lo) are real 
at ujn — 0. Comparing Eq. (|26|) with the ujn — 
term of Eq. H48I) . we see that the kernel coefficients re- 
duce to the time-independent perturbation theory ex- 
pressions for the second derivatives of the total energy 
with respect to exchange field, d'^Etot{^MF)/dA^ dA^ , 
i9^£^tot(AMF)/A^ dAy. In its static limit, the fluctuation 
kernel reduces to one that would be obtained from a clas- 
sical theory with the micromagnetic energy functional 
derived from a mean-field-theory calculation, or for accu- 
rate first principles calculations by solving spin-density- 
functional Kohn Sham equations self-consistently. By ex- 
panding Kapiiujn) around iuJn — 0, we obtain 
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Kyy{iUJn) 



MF) 



9'i;tot(A 



MF) 



dl^y dAv 



d^EtotiA 



MF) 



dA^ dAv 

d"^ Etot{AyiF) 

dAv dAv 



+ Oiiul)^ia + b) + Oiiul), 
+ OiLul)^{a~b) + OiLul), 
- C(Amf) iujn^ c- C(Amf) iuj„ 

+ C(Amf) it^n = C + C(Amf) ioJn 
I 



(60) 
(61) 
(62) 
(63) 



where 



U 2 



are real constants and the Berry curvature is 



C(Amf) = ^E 



2^ (6-6)^ 



(64) 



(65) 



(66) 



(67) 



This expression for the Berry curvature can also be 
derived from Eq.( I29|l by using time- independent per- 
turbation theory expressions for the dependence of 
single-particle wavefunctions on exchange-field orienta- 
tion. Note that the leading frequency dependences in 
Kxx and Kyy are quadratic. These equations are valid 
for iujn smaller than the smallest particle-hole excitation 
energy — as we discuss at greater length below. 
So far we have completely disregarded amplitude fluctu- 
ations of the exchange field,?. e. the components of the 
Kernel involving A^. It is obvious that in absence of 
spin-orbit interactions 



Ky z{iuJn) 



KzxiiUJn) = , 
Kzy{iuJn) = , 

A/a Snfl_ sr^ dn-p 
U 4 ^ 



(68) 
(69) 

(70) 



When spin-orbit is present, the off-diagonal components 
involving are in principle non-zero. However it turns 
out that they are always very small (see next section); 
therefore we will typically neglect them and we will only 
keep Kz z{i'^n) = J^a/U. This approximation amounts 
to neglecting the weak dependence of the magnitude of 
the exchange field on its orientation. 



C. Excitation Spectrum 

The kernel of the fluctuation action is the matrix in- 
verse of the exchange-field propagator, which is propor- 
tional to the spin susceptibility linear response func- 
tion. We refer to this propagator below as the spin- 
susceptibility. 



(71) 



The elementary magnetic excitations of the nanoparti- 
cle (particle-hole and collective) occur at real frequencies 
frequencies uj at which Xa/3(w) has poles. Here Xa/3(w) is 
obtained by analytically continuing Eq. (|71f) to real fre- 
quencies: iLOn Lo + irj. Let D{u}) denote the complex 
function, 



D{u;) = det[K{Lj)] 



(72) 



The condition for the existence of a pole in Xafiij-^) can 
be expressed by the equation 



D{uj) = 



(73) 



In order to make progress, we will assume that all the 
components of K involving the longitudinal variable z 
are small and can be disregarded, except for zi^) The 
expression for the elementary excitation energy then sim- 
plifies to 

D{lj) w Kzz[Kxx{cj) Kyy{u) - Kxy{uj) KyxiLo)] = 0(,74) 

where we have used Eqs. (|53l I55|) . 

Let us consider first the limit of small nanoparticles, 
where the single-particle mean-level spacing 5 is much 
larger than the ferromagnetic resonance, which is of order 
the anisotropy energy per atom, K. In Cobalt K k, Q.l 
meV, and if use bulk single-particle density of states [l7| 
we find that this limit is reached in nanoparticles con- 
taining Na « lO'' atoms. In this case we expect a 
pure ferromagnetic resonance mode with a large spectral 
weight that appears as a separate quantum state below 
the lowest particle-hole excitation energy. (We discuss 
the situation where particle-hole and collective excita- 
tions are not cleanly separated at greater length below.) 
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we obtain 

0, (75) 



Using the expansions of Eqs 

Ma 

U 

which yields a low-energy pole at the real frequency 
^/ a? ~ b"^ — 



Around the pole we have, 



|iC(AMF)| 



(76) 



-1 -Zio 



where the "residue" Z^^s is 



= 2 ^ Va^ - 62 - 



UC(A 



(77) 



(78) 



Thus the collective excitations have a gap at E'rcs- The 
gap is proportional to the quotient of ^ a?' — h"^ — , 
which is essentially the total anisotropy energy, and the 
Berry curvature C(Amf)- For particles that are not too 
small, both quantities are proportional to the particle 
volume, thus i?res is approximately independent of par- 
ticle volume and approximately equal to the anisotropy 
energy/atom, K. The spectral weight of the collective 
mode Z-cca is inversely proportional to the product of the 
same two quantities. The spectral weight divided by the 
resonance frequency is proportional to the static response 
of the exchange- field, a quantity which can be understood 
simply by minimizing the total micromagnetic energy. 

The considerations described above apply when the 
collective excitation energy is smaller than the lowest en- 
ergy particle-hole excitation. Even in this limit where 
particle-hole and collective excitations are well separated, 
the exchange-field propagator does have some spectral 
weight in particle hole excitations as well. The Kernel 
(Eq. I49|l must have zeroes between the poles that occur 
at each mean-field particle-hole excitation energy. There- 
fore we expect D{uj) = at w = Coji « uoji = — Cj ^-nd 
the susceptibility to have additional poles at cjji. Let us 
expand D{uj) around one of these poles 



(79) 



Thus, in general the residue of each pole is proportional 
to 



(80) 



D. 



Coupling between Collective and Particle-hole 
excitations 



If we increase the size of the nanoparticle to a few 
thousands atoms, the lowest particle-hole excitation en- 



In this situation the pole at E'res will start to loose its 
collective character. As 5 becomes much smaller than 
ii^rcs, the ferromagnetic resonance mode will no longer 
be a single excitation, and instead appear as enhanced 
spectral weight in response functions that is spread over 
several particle-hole excitations. If we further increase 
the nanoparticle size such that 5 << -Eres, we start to 
approach the thermodynamic limit, where the collective 
ferromagnetic resonance will be spread over a large num- 
ber of particle-hole excitations. If i^ics is the frequency 
at which Eg. 1741 holds, by expanding D{ljj) near Ey-^s the 
ferromagnetic collective mode has a width given by 



r 



, Im D{E,,,) 
'ReD'{E,,,) 



(81) 



The susceptibility then has a resonant denominator and 
ferromagnetic mode has the shape of an asymmetric 
LorentzianpGj. 

It is possible to relate the typical value of matrix el- 
ements between single-particle states that appear in the 
kernel for a large nanoparticle and the Gilbert damp- 
ing parameter usually used to characterize the width 
of a ferromagnetic resonance line. The Gilbert damp- 
ing parameter a is the ratio between the line-width 
and the resonance frequency. It is normally introduced 
as a phenomenological damping parameter in Landau- 
Liftshitz equations of motion for the magnetization di- 
rection; these equations are implied by the low-frequency 
dynamics discussed above. Taking the continuum limit 
for the particle-hole excitation spectrum and assuming 
that we have approximate rotational symmetry about 
the exchange-field-orientation axis, it follows from the 
preceding analysis that [s^ 



[UJ 



E„ 



a 



llR[K,y{LU = Sres)] 252|C(Amf)| 



\{J\-s.\lW 



ergy ujj 



5 will start to approach E'l-cs from above. 



(82) 

where the overbar denotes an average over typical 
particle-hole matrix elements for states near the Fermi 
energy. Note that a approaches a constant and the typ- 
ical matrix element scales like ^J'X^ in the limit of large 
nanoparticles. 

We can estimate the number of individual particle- 
hole excitations that contribute to the ferromagnetic res- 
onance in a large nanoparticle. The density-of-states 
for particle-hole excitations grows linearly with energy 
and at the resonance energy is ~ E,-cs/S'^. The num- 
ber of particle-hole excitations in the resonance is the 
width of the resonance aE^cs times this density of states, 
~ aE'^^^/5'^. This expression implies a condition for 
the crossover to pure-state ferromagnetic resonance when 
5 — \foL Ercs- The crossover is therefore expected to oc- 
cur at larger particle sizes when the ferromagnetic reso- 
nance is sharp. Even when the particle is small enough 
that the resonance is not normally coupled with particle- 
hole excitations, avoided crossings between these two 
types of excitations will frequently occur as the exter- 
nal magnetic field or other parameters are varied. We 
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can describe such an avoided crossing by assuming that 
the low-frequency Hmit can be taken for all but one of 
the particle-hole excitation contributions to the kernel. 
To briefly explain what happens in this limit, for which 
the kernel can be written as the sum of a part if™""*'' 
with a smooth frequency dependence and a part K^^°^ 
with a resonant frequency dependence, we assume that 
the particular particle-hole excitation contributes only to 
Kl°^, making a contribution 

K^:,%u;)^^^\{J\sM)\'- (83) 

LU — LU ■ ■ 

It follows that the poles of the exchange-field propagator 
occur at energies 

^± = ^^^^±^ ± [m^. - u,,)/2Y + V'] (84) 

where V is an avoided crossing gap. Using Ea.f I82|) we 
find that V ^ aE'rcs- The size of the avoided crossing gap 
that occurs in a small nanoparticlc when a particle-hole 
excitation energy is tuned through a collective excitation 
energy by an external field or another parameter is spec- 
ified by the ferromagnetic resonance width in the limit of 
a large particle. 

It is important to note that the experiments in Ref.[3,0 
are carried out for nanoparticles containing Ma ~ 1500 
atoms. Thus these experiments are probing the most 
interesting and difficult intermediate regime of parti- 
cle size, where the lowest particle-hole excitation energy 

V. NUMERICAL RESULTS AND DISCUSSION 

In this session we will present numerical calculations 
performed on our two ferromagnetic nanoparticle models 
which illustrate and support all the main points of the 
theory of elementary spin excitations developed above. 
Since we can deal with nanoparticles containing up to 
260 atoms, we expect, on the basis of our theoretical con- 
siderations, that the exchange-field correlation function 
will normally display one main peak at energies below 
the lowest particle-hole excitation. In addition to verify- 
ing this expectation, we will see that the ferromagnetic 
resonance mode is characterized by a very interesting be- 
havior when manipulated with an external magnetic field 
to bring it through con incidence with a particle-hole ex- 
citation. 

As mentioned already in Sec. II VI in discussing the dif- 
ferent components of the kernel and the susceptibility we 
will take the z-direction along the direction of the mag- 
netization, fi. We will also choose the y-dircction to lie in 
the equatorial XY-plane of the nanoparticle. Then the 
x-direction will be fixed by the condition of being orthog- 
onal to both z and y. Therefore the x and y components 
of the kernel describe transverse spin fluctuations with 
respect to the direction of the magnetization, but these 




FIG. 2: Anisotropy landscape i5tot(A) as a function of and 
$ for a 143-atom hemispherical nanoparticle. -Etot is periodic 
in <I>, with period n/2. 



directions in general will not be symmetry directions with 
respect the crystal structure or the nanoparticle geome- 
try. The direction of Cl is determined by minimizing the 
the classical micromagnetic energy functional -Etot (A), 
or equivalently by solving the mean-field equations self- 
consistently. The variation of Etot (A) as a function of 
& and $ has been studied in detail in Ref. |3|. The 
main conclusions are summarized in Fig. |21 where we 
plot Etot (A) in the (0, $)-plane at /?oxt = for a hemi- 
spherical nanoparticle consisting of Ma = 143 Co atoms, 
arranged in a f.c.c. lattice. In this case the magnetization 
direction f2 lies in the XY-plane (0 = 7r/2), along one of 
the four directions corresponding to the four degenerate 
shallow minima of Etot (A). By applying an external field 
some of these minima will become classically metastable 
and we can have different hysteretic behaviors with zero, 
one or two reversals of the magnetic moment, depending 
on the direction of the field. 



A. Real part of Kxx(io) 

We first consider the spectral representation of the ker- 
nel defined in Eq. In Fig.|3Ia) we plot ReK^xii^^n) 
vs. Re(«(x'„) ~ after analytical continuation iun — > uj+ir] - 
for a 26-atom nanoparticle for the case of the LRE model. 
Consistent with Eqs. l(S7l IS^ . for zero SO we find that 
the Kernel has just one pole at particle-hole excitation 
energy equal Amf, and a gapless zero. With SO interac- 
tion included, many more matrix elements (J|s^|/) will 
be non-zero, and consequently many more poles corre- 
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FIG. 3; Real part of the transverse diagonal component of the 
kernel for a hemispherical 26-atom nanoparticle, Long-Range 
Exchange model, (a) Without spin-orbit (SO) interaction, 
(b) With SO interaction. 



spending to these particle-hole excitations will appear in 
the kernel. This is shown in Fig. |2Ib), where we can 
see that the poles with largest residues still correspond 
to particle- hole excitations near Amf- Although not vis- 
ible in Fig. Olb) because of its very small residue, the 
first pole in the kernel occurs at an energy of order of 
the single-particle mean level spacing (5, which for this 
particle size is approximately 20 meV. This is shown in 
Fig. ^ where we zoom on a small energy window at very 
low energies. 

Similar results for the LDE model are shown in 
Fig. (0). We can see [Fig. El^a)] that the kernel has al- 
ready several poles even in absence of SO coupling, since 
not all the energies of majority- and minority-spin states 
are shifted by the same rigid spin-splitting field, as they 
were in the case of the LRE model. Most of the poles 
are still concentrated near Amf, however. The special 
properties of the LRE model which result in rigidly split 
mean- field bands do not, therefore, introduce and special 
artificial features in the low-energy excitation spectrum. 
From Fig. Efb) it would seem that SO interaction do 
not introduce major effects, merely broadening the pre- 
existing pole structure. In fact SO interaction changes 
the analytical structure of the kernel at low energies, as 
it did for the LRE model. Indeed, Fig. shows that 
when SO is included the first pole of K^x occurs again 
at an energy of the order of the single-particle mean level 
spacing. We find that most of these low energy particle- 
hole excitations involve pairs of states that have mostly 
the same minority spin character, as expected because 
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FIG. 4: Low-energy behavior of the Kernel for the same sys- 
tem as in Fig. |3| Solid line: without SO interaction; dashed 
line: with SO interaction. The vertical dotted line marks the 
position of the first particle-hole excitation, where YieKx x has 
a pole. 
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FIG. 5: Real part of the transverse diagonal component of 
the kernel for a hemispherical 26-atom nanoparticle. Local 
D-orbital Exchange model, (a) Without SO interaction, (b) 
With SO interaction. 
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FIG. 6: Low-energy behavior of the Kernel for the same sys- 
tem as in Fig. |K| Solid line: without SO interaction; dashed 
line: with SO interaction. The vertical dotted lines mark the 
position of particle-hole excitations. 



FIG. 7: Real part of the transverse diagonal component of the 
kernel at low-energies, for a hemispherical 143-atom nanopar- 
ticle. (a) Long-Range Exchange model, (b) Local D-orbital 
Exchange model. Vertical dashed lines mark the position of 
particle-hole excitations, where Reifx ^ has poles. 



essentially all states around Fermi level are minority- 
spin states. Obviously matrix elements ( J, | |s+|/, |) 
are identically zero and therefore there are no poles at 
these low energies when SO coupling is not present. It 
is only when SO coupling is present that these states ob- 
tain a small admixture of majority-spin character which 
produces non-zero matrix elements. It is important to 
note that even when the rare presence of majority spin 
states near the Fermi level is recognized, their matrix 
elements with minority spin states close in energy are al- 
ways very small because their orbital wavefunctions are 
almost orthogonal. Fig. [7| shows the low-frequency be- 
havior of the kernel for our two models for the case of a 
143-atom nanoparticle. We see essentially the same trend 
as for the case of smaller nanoparticles: when SO is in- 
cluded the first pole of the Kernel occurs at an energy of 
the order of the single-particle mean-level spacing, which 
for this nanoparticle size is of the order of a few meV's. 
Note that also at this particle size the residue of this low- 
energy pole is three orders of magnitude smaller than the 
residues of the poles occurring near Amf- 

The existence of spin-flip particle-hole excitations at 
energies of order b has important implications: we will 
show below that these particle-hole pairs can couple to 
the low-energy spin collective mode when, by increas- 
ing the nanoparticle size and/or by applying an external 
field, their energy start to approach i?rcs- 



B. Spectral function: Imxia;(ii') vs. u) 

We now discuss the spin-fluctuation spectral function, 
that is the imaginary part of the spin-susceptibility de- 
fined in Eq. ifTT)) . In Fig.|Slwe plot \m.Xxx{!^) as a func- 
tion of a; for a 143-atom hemispherical nanoparticle for 
several values of the external magnetic field. As expected 
form our discussion in Sec. IIVI for this nanoparticle size 
we find that in general the spectral function has only one 
pole that has substantial weight. At low external fields 
this pole occurs at a frequency that satisfies Eq. H74|l . well 
below the lowest particle-hole excitation. It can therefore 
be identified as a spin-collective mode. In this regime 
the collective mode is an exact elementary excitation, at 
least within the Gaussian approximation that we con- 
sider. The finite width in Fig. |Slis due to a finite 77 in 

iiOn ^ Ul + it]. 

Apart from the ferromagnetic resonance pole, the sus- 
ceptibility has other poles at higher energies, as seen in 
Fig. El where we plot D{uj) given in Eq. These poles 
are all very close to MF particle-hole excitation energies 
and have virtually zero spectral weight for this particle 
size, when the external field is zero. We can therefore 
identify them as the spin-flip particle-hole excitations dis- 
cussed in Sec. IIVCI 

As shown in Fig. |S1 we can use an external magnetic 
field to manipulate both energy and spectral weight of 
the collective mode. From zero field to the reversal field, 
the collective mode energy and its spectral weight de- 
crease. Beyond the coercive field they both start to in- 
crease monotonically. As shown in Fig. |Htb) , if we keep 
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FIG. 8: Imaginary part of the transverse diagonal spin sus- 
ceptibility {x component) at different external magnetic fields, 
for a hemispherical 143-atom nanoparticle in the LDE model. 
Hext is in the ZX -plane at an angle 7r/4 with the Z-axis. 
J^ext = 0.15 T is a reversal point. This low-energy peak in 
the susceptibility corresponds to a spin-collective mode due 
to coherent magnetization fluctuations. Vertical lines in (b) 
mark the position of the first particle-hole excitation energy, 
which shifts down with increasing field. When the collective 
mode energy approaches the particle-hole excitation energy, 
the collective mode peak splits and part of its spectral weight 
is taken away from the pole. 



increasing the field, the collective mode energy starts to 
approach the first particle-hole excitation energy, whose 
energy decreases with increasing field. When the two en- 
ergies are proximate to each other, the collective-mode 
peak splits into two close resonances: part of spectral 
weight of the magnetization fluctuation is transferred to 
the nearby particle-hole excitation, as also illustrated in 
FigllOl where we plot the derivative of D{w) at Ej-cs as 
a function of the external field. The peak in D' corre- 
sponds to a sudden decrease of Zj-cs given in Eq. IHOl the 
collective mode can decay into a particle-hole excitation. 
More precisely, we can say that the two types of exci- 
tations are coupled, and therefore it is no longer mean- 
ingful to speak about spin-collective modes as distinct 
from spin-flip particle-hole excitations. In our model the 
coupling mechanism is provided by the spin-orbit inter- 
action. Indeed, we have checked that without spin-orbit 
interaction, nothing happens to the collective peak when 
its energy crosses a particle-hole excitation. 
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FIG. 9: Determinant of the kernel defined in Eq. 1741 for 
the system of Fig. |H| at zero external field. The first zero of 
KeD{u>) at 0.06 meV corresponds to the ferromagnetic res- 
onance seen in Fig. |H] KeD{u}) has other zero's very close 
to its poles that are located at particle-hole energies (marked 
by vertical dashed lines). The zero's correspond to spin-flip 
particle-hole excitations, but their spectral weight in Fig. |H| 
is virtually zero unless they approach the ferromagnetic reso- 
nance mode. 
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FIG. 10: Derivative of the determinant of the kernel given in 
Eg. 1741 calculated at Eros, as a function of the external field for 
the system of Fig|H| D'{Ercs) is inversely proportional to the 
residue of the collective mode. The peak in D'{Eies) occurs 
where the collective mode energy and the first particle-hole 
excitation energy cross. 
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FIG. 11: Variation of the collective-mode energy (low-energy 
peak in the susceptibility) as a function of the external field 
for 143-atom and 260-atom nanoparticles. //ext is in the 
ZX-plane at an angle 7r/4 with the Z-axis. In (b) the de- 
tailed low-field behavior is shown. The discontinuities for the 
143-atom nanoparticle correspond to magnetization reversal 
points. (There are no reversal points in the 260-atom system 
for this external field direction.) 



The coupling between collective modes and particle- 
hole excitations described above mimics vifhat would hap- 
pen if, instead of manipulating the excitation energies 
with an external magnetic field, we could progressively 
increase the size of the nanoparticle up to a few thou- 
sand atoms. In this case -Eres would stay approximately 
constant but S would decrease. When these two energy 
scales cross the lowest particle-hole excitations will start 
to interact with the collective mode. It is therefore ex- 
pected that in a nanoparticle with a few thousand atoms 
where i?ros ~ 5, spin collective modes and particle-hole 
pairs are strongly coupled. 

Fig. II II shows the external field dependence of the en- 
ergy of the collective mode peak for a 143-atom and a 
260-atom nanoparticle respectively. The jumps in the 
energy for the 143-atom nanoparticle correspond to re- 
versal of the magnetic moment. Interestingly, we have 
found that the peak energy follows accurately Eq. [THI for 
all values of the external field, even when its value is 
larger than the lowest particle- hole energy [s^. The en- 
ergy gap Ej-es, obtained extrapolating the curve at zero 
field, is approximately volume independent as expected 
from our general theory. We would like to point out that 
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FIG. 12: Imaginary part of the transverse diagonal spin sus- 
ceptibility, X and y components. The external magnetic field 
is in the ZX-plane at an angle 7r/4 with the Z-axis. At 
high external fields the two transverse components of the 
spin susceptibility become equal, (a) Hemispherical 143-atom 
nanoparticle. (b) Hemispherical 260-atom nanoparticle. 



at low fields (i?cff < IT) i^rcs ~ 0.1 meV, which of the or- 
der of the anisotropy constant /atom K in Cobalt 8] , and 
is also of the order of the tunneling resonance spacing 
observed in Refsllll. 



C. Anisotropic Fluctuations 

The dynamical susceptibility displays spatial 
anisotropy in its transverse diagonal components 
due to the strong anisotropy of the micromagnetic 
energy functional. We illustrate this point in Fig. E| 
by plotting hnxxxi^) together with Im Xyy{(^)- We 
can see that also Imxyyi^^) has only one dominating 
collective-mode peak at the same energy of the collective 
mode in Imxxxi^)- This is obvious since in both cases 
the collective mode energy is given by Eq. H74|) . However 
the weight of the pole in the two spectral functions is 
different and its variation as a function of the external 
field opposite. At large fields the difference between 
X and y components disappears and the transverse 
susceptibility becomes isotropic. We find the same trend 
for hemispherical nanoparticles containing 143 or 260 
atoms. 

We can get an intuitive understanding of this behav- 
ior by looking again at the classical micromagnetic en- 
ergy functional Etot{^) as a function of 8 and (f>, shown 
in Fig. refanisotropy. There are four equivalent min- 
ima in the XY plane at zkX ± Y, that is = 7r/2, 
$ = 7r/4 + n7r/2, n = 1,2,3. Note however that the 
energy barrier separating the four minima at Q — tt/2 
is very low. We can interpret the collective mode as 
the zero point motion of a two-dimensional anisotropic 
harmonic oscillator, whose potential is obtained by ex- 
panding i?tot(0,*l') to second-order in 8 and $ around 
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one of the minima. Looking at the x-component of the 
spectral function corresponds to exciting the collective 
oscillation mainly in the hard direction. We can pursue 
further this analogy and clarify the effect of the external 
magnetic field. A magnetic field in the ZX plane starts 
to decrease the low barriers separating the four minima. 
Thus the collective mode energy decreases. This con- 
tinue up to reversal. Note that at reversal the collective 
mode energy does not go to zero because the potential 
landscape has only a saddle point there: the potential 
still increases in the direction of the two poles of the 
unit sphere {Q = 0, 7r/2). After the last reversal has 
taken place, namely once the magnetization has flipped 
to a stable minimum, a further increase of the magnetic 
field makes the spring constant of the harmonic oscilla- 
tor stiffer. Hence the collective mode energy starts to in- 
crease. This behavior is summarized in Fig llll At strong 
magnetic fields the energy functional is dominated by the 
Zeeman term and any anisotropy in the transverse fluc- 
tuations disappears. 

Spherical particles have a more symmetric anisotropy 
energy landscape, without the strong XY easy plane 
anisotropy characteristic of hemispherical particles. A 
reversal point in this case will correspond the disappear- 
ance of the energy barrier in £^tot(^^) in all directions per- 
pendicular to Amf, or at least a less pronounced saddle 
point. Therefore spherical nanoparticles show a much 
stronger decrease of collective mode energies at reversal 
points. 



D. Berry Curvature 

We conclude this session with a discussion of the nu- 
merical results for the Berry curvature C(Amf) defined 
in Eq. We have seen that C(Amf) affects the quan- 
tization condition of the collective mode energy given in 
Eq.[7ni and is inversely proportional to its residue. Eg 1781 
Variations of C(Amf) from the constant value iS/A^^p 
reflect the non-trivial role played by the spin-orbit inter- 
action. In Fig. El we plot C(Amf) as a function of an 
external magnetic field for a 143-atom nanoparticle. We 
have computed C(Amf) in two different ways, according 
to Eg.lHTIand Eas . I34I3'51 respectively. The agreement be- 
tween the two calculation methods is excellent, especially 
at low fields [see fig. I13f b) ] . The fact that the compu- 
tation method based on Eq. 1351 works well is significant 
because this method relies only on the knowledge of the 
ground-state for a given magnetization orientation, which 
in principle can be obtained with fairly good accuracy 
from density-functional calculations. Thus Eq. 1351 pro- 
vides a very convenient expression for computing C(Amf) 
beyond the mean-field approximation used in the present 
work, and also for more realistic nanoparticle models. 

From Fig. El we can see that the external field depen- 
dence of C(Amf) is rather smooth and weak, except when 
the system approaches a reversal point. At a reversal 
point C(Amf) suffers a discontinuous jump. On the other 




FIG. 13: Berry curvature C(Amf) for a 143-atom nanoparti- 
cle, computed using two different methods: the solid line is 
obtained from Eq. 1671 the dashed line from Eqs. I34I35I (b) 
Behavior at low fields. 



hand, C(Amf) is completely insensitive to the crossings 
mentioned above between the ferromagnetic resonance 
mode and particle-hole excitations, when the latter are 
of order « S. In fact, if we look at the perturbative 
expression of C(Amf) given in Eq. 1671 we can see that 
there is no reason to expect any large fluctuations of this 
quantity when the smallest energy denominator in the 
sum is « S. If the collective mode energy approaches a 
particle-hole excitation, C(Amf) will still be well-defined 
and smooth, but it will not tell us anything about how 
the collective mode spectral weight is depleted in favor 
of the nearby particle-hole excitation. In other words, 
the residue of the collective mode is inversely propor- 
tional to C(Amf) only when the low- frequency expan- 
sion is valid, in which case it will be a smooth quantity. 
When the collective mode energy is close to a particle- 
hole excitation energy of order 6, the low-frequency ex- 
pansion is meaningless. One must look at the derivative 
of the determinant D{u;) calculated at the pole frequency, 
given in Eq. I8UI to find the residue. An important ex- 
ception occurs when a particle-hole excitation energy ap- 
proaches zero (i.e. close to mean-field quasiparticle level 
crossings) 14]. These are events that can occur only for 
some particular orientations of the magnetization and at 
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isolated values of the external fields. When the system 
is close to a mean-field quasiparticle level crossing, the 
Berry curvature landscape C(A) is strongly distorted. In 
this case the fluctuations of C(A) are a direct indication 
of the coupling between the collective mode and particle- 
hole excitations, which obliterate the distinction between 
them. 



VI. CONCLUSIONS 

In summary, we have developed a theory of elemen- 
tary spin excitations in ferromagnetic metal nanopar- 
ticles that provides a consistent and unified quantum 
description of both quasi-particle and collective mode 
physics. Our formalism, based on a path integral ap- 
proach, allows us to make a connection between micro- 
scopic exchange and spin-orbit interactions and classi- 
cal micromagnetic theory. We have shown that small 
nanoparticles have a collective excitation at energies be- 
low the lowest particle-hole excitation energy, whose en- 
ergy gap i?ros is given by the ratio of the anisotropy en- 
ergy and the total Berry phase of the system. As the 
single-particle mean-level spacing decreases and becomes 
much smaller than ^/aE^as, the collective excited states 
evolves into a damped collective mode whose spectral 
weight is distributed over a large number of particle-hole 
excitations. 

We have illustrated these ideas by performing numer- 
ical calculations of nanoparticles containing up to 260 
atoms, described by a microscopic tight-binding model. 
We have found that for this particle size there is typi- 
cally an isolated collective excited state below the lowest 
particle-hole excitation which nearly exhausts the spec- 
tral weight of the dynamical susceptibility. The energy 
gap Ej-Qs is of the order of O.lmeV, and is approximately 
independent of the particle size. Occasional crossings be- 
tween primarily single-particle and collective excitations 
occur as a function of applied magnetic. Near the cross- 
ing point, the collective mode peak splits, as a result of 
resonant coupling between the two types of excitations 
that is non-zero because of spin-orbit interactions. These 
crossings become more common as the particles become 
larger and cannot be avoided for system for nanoparticles 
with more than typically 10,000 atoms. 

Although a detailed comparison with the tunneling 
transport experiments of Ref. 0, is beyond the scope 
of the present paper, our analysis sheds light on some 
essential features which are quite relevant to the under- 
standing of the experimental results. What emerges from 
our theory is a picture of elementary excitations that is 
far more complex than that derived from earlier phe- 
nomenological models, where the effect of spin-orbit in- 
teraction is accounted for only indirectly by means of 
an uniaxial anisotropy term in a giant spin Hamiltonian 
that represents the coherent magnetization dynamics. In 
particular, our results suggest that for the particle size 
considered in the experiments quasi-particle and spin col- 



lective modes are most likely strongly entangled by spin- 
orbit interactions. It is not unconceivable that such an 
intertwined set of excitations could provide a rich low- 
energy tunneling spectrum, even in conditions of equilib- 
rium tunneling. 

The energy gap of the ferromagnetic resonance mode 
can be viewed as the characteristic mean-energy-level 
spacing between coupled collective-quasiparticle excita- 
tions described by an effective Hamiltonian. Interest- 
ingly enough, the value of i^rcs deduced from our model 
calculations is of the order of the observed tunneling reso- 
nance level spacing. This value, which is also of the order 
of -ftTbuik, is five times larger than the anisotropy con- 
stant estimated from the measured switching field using 
Ksw = Hb -ffsw • The smallness of Kgw is one of the prop- 
erties that led the authors of Ref@ to conclude that equi- 
librium spin excitations involving only the lowest spin- 
multiplet cannot be resolved in the present experiments 
and therefore cannot explain the observed large density 
of resonances. While the discrepancy between Kgv, and 
-K'buik is still a puzzle (see however [sj), our results 
suggest that the these low-lying excitations are in fact be- 
ing detected. These considerations do not necessarily im- 
ply that the main features of the tunneling experiments 
can be understood by equilibrium transitions alone. Mea- 
surements on gated devices 3] convincingly support the 
hypothesis that non equilibrium transitions play a crucial 
role. The point that we want to make here is that low- 
lying equilibrium spin excitations are probably equally 
important for the interpretation of the experimental re- 
sults. 

It would be highly desirable to perform new tunnel- 
ing experiments in more controlled situations and espe- 
cially with smaller nanoparticles, so that the spin exci- 
tations associated with the electronic degrees of freedom 
and with the magnetization collective coordinate could 
be more easily disentangled. The physical condition re- 
quired to reach this regime is that the single-particle 
mean level spacing be larger than the total anisotropy 
energy jl3|. For Cobalt this would imply dealing with 
nanoparticles containing of the order of 100 atoms, which 
is within reach of the present experimental capabilities. 
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